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Figure Captions 



1. Picture of a cluster grown using conformal growth model. 

2. Illustration of effect of / on the unit circle. 

3. Diagrams illustrating interaction vertex and contraction of noise terms. The first diagram 
results from setting G{j) = W[j, G{k), G{1), G* (m)]. The second diagram would occur in evalu- 
ating {Gf{j)Gf*{k)) = 5{j - k)NU). 

4. a) A possible contribution to the average of b)A possible contribution to s'^-^^. c)A 
diagram which should not be considered if the propagators are already drawn in terms of s^-^-^. 

5. Perturbation expansion in the example theory. 

6. Diagram contributing to the renormalization of s^-^-^. 

7. a) Diagram contributing to the renormalization of t. b)Another possible contribution to 
the renormalization of t. This contribution is ignored. Certain difficulties are encountered in 
drawing this diagram. If the t interaction were drawn in standard form, many lines would need 
to cross each other. 

8. Other possible contributions in the RG flow. 

9. a)Simplest diagram for the correlation function needed to compute t(5). The external lines 
are the various G{j), G*{k) in the correlation function, b) Tree diagram for the same correlation 
function. 

10. Plot of scaling of A against cutoff, as described in text. 

11. Plot of scaling of log of mean square of G{j) against log of j, as described in text. 
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An analytical renormalization group treatment is presented of a model which, for one 
value of parameters, is equivalent to diffusion limited aggregation. The fractal dimension of 
DLA is computed to be 2 — 1/2 + 1/5 = 1.7. Higher multifractal exponents are also calculated 



I and found in agreement with numerical results. It may be possible to use this technique to 

. describe the dielectric breakdown model as well, which is given by different parameter values. 



1 Introduction and Notation 



C , 1.1 Introduction 

o 



Diffusion limited aggregation (DLA) is a model for growth of a cluster Q, by the accretion of 
random walkers. These random walkers arrive from infinity and stick when they contact the 
cluster. After a walker sticks to the cluster, the next walker is released from infinity. This 
process gives rise to fractal patterns. Due to the mathematical equivalence of random walks and 
potential theory, this procedure is equivalent to solving Laplace's equation outside the boundary 
of the aggregate, setting the potential zero on the aggregate and constant at infinity, and picking 
a point on the surface of the aggregate to add the walker with a probability proportional to the 
local field strength; this field strength may be thought of as an "electric field" . 

There has been much numerical work on DLA in two dimensions, where the fractal 
dimension has been determined to be 1.71^]. Mean-field calculations predict D = 5/3^, 
which indicates that, in two-dimensions, something is lacking in the mean-field approach. In 
higher dimensions the mean-field theory appears much more accurate. One first principles 
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renormalization theory, based on the branching nature of DLA processes, obtained the result 
that D = 1.661 @. 

Another important analytic result was the derivation of the electrostatic scaling law, 
which appears to be obeyed numerically by the aggregates [^. This law is used as an essential 
step in the calculation of this paper. 

Recently, another formulation of DLA was proposed, based completely on conformal 
mappings[^. A conformal function mapped the unit circle onto the boundary of the aggregate. 
In this formulation of the problem, the electrostatic scaling law followed almost automatically 
when considering the behavior of the first Fourier component of the mapping. We investigated 
numerically the problem of the importance of different Fourier components of the mapping. It 
appeared that by directly simulating the dynamics of only on the first few Fourier components, 
results could be obtained for short growth periods that were similar to those when the full 
function was used. This suggested that it might be possible to develop a renormalization theory 
based on integrating out higher Fourier components, using techniques similar to those used in 
field theory. 

Also, by comparing a picture of the cluster generated by only keeping some small number 
of terms in the Fourier expansion of the mapping to a picture generated by the full mapping, it 
appeared that the finite number of terms gave a good description of the boundary of the object. 
It did not accurately describe the exact microscopic structure of the growth tips, and did not 
correctly describe the structure of portions of the object far from the growth region, that is, 
deep in the inside of the object. However, one would expect that the microscopic structure is 
not too important, and that the description of regions where there is little probability of growth 
is also not important. 

Fig. 1 shows a picture of the cluster that results from the conformal mapping model. 
The envelope surrounding the cluster was generated from the first 40 terms of the Fourier series 
expansion of the mapping used to generate the full cluster. 

For longer time periods, more terms in the Fourier series were needed, but this is only 
to be expected; if only one term were kept in the Fourier series the object would be a circle and 
would grow with a radius proportional to the square root of time. As more terms are kept, the 
object can grow faster than the square root of time by changing shape, but for any given number 
of terms, eventually the growth will be as the square root of time. Therefore, it is expected that 
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in the RG that follows there will be some cutoff in the number of terms kept which increases 
with the size of the object. 

The above discussion is intended to motivate the RG that follows. Most of the discussion 
is done in more detail in Ref. Q. 

Hope, that such a scheme would work, was provided by numerical evidence that the 
conserved quantities (the moments) of the continuum growth law were very nearly conserved 
by the random growth process^. For the lowest Fourier coefficients, one would expect that the 
random growth would be close to the average growth determined by the continuum law, while 
the higher Fourier coefficients would fluctuate more wildly. 

The paper is divided into several parts. First, the conformal mapping model for DLA 
is discussed and used to derive continuum equations for growth, essentially equivalent to the 
Shraiman-Bensimon equation for the Hele-Shaw problem These equations are heuristically 
modified to add the essential differences between DLA and its continuum limit: the presence of 
noise and the existence of a microscopic cutoff. This leads to a new model which is hoped to 
be in the same universality class as DLA. Even if it is not in the same universality class, it is 
similar enough to be of interest in itself. 

Second, under an adiabatic assumption, the equation is modified to vastly simplify the 
time dependence, leaving almost a static problem. The adiabatic assumption makes possible the 
RG and perturbation theory calculations described latter in the paper. The adiabatic assumption 
is justified by numerics and self-consistently by the RG itself. 

At this point, before doing the RG, it is still possible to make some comparison to 
numerics based on the continuum representation of DLA. 

Third, a perturbation theory is developed for the continuum equation, with a well de- 
fined set of rules for calculating correlation functions. The perturbation theory requires some 
resummation of diagrams, where to calculate resummed propagators it is necessary to use a 
renormalization group approach. This RG forms the fourth part of the paper; the calculations 
for the RG have only been done to lowest order, producing an appropriate renormalized propa- 
gator and vertex. Fifth, the adiabatic assumption is removed, and the renormalized propagator 
is used to calculate various exponents in the theory. Sixth, the results are compared to numerical 
experiments, the self-consistency of the approach is discussed, and there is discussion of what 
may happen if the computation is performed to higher orders. 
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1.2 Notation 

A large number of functions will need to be defined in this paper. As much as possible, I will 
use the following notations. Capital letters are used for functions, such as F, G to be defined 
latter, which describe the shape of a specific growing aggregate. Power series expansions of 
these functions will be denoted by subscripts, so F{z) = Fiz^ + Fqz^ + F_iz~^ + ... In the 
continuum limit of these power series expansions, to be appropriately defined latter, where sums 
are approximated by integrals, the Latin letters j, fc, /, m, n, o will be used as indices. One will 
see terms like G{j). 

Greek letters e, A, /i will be used for ultraviolet and infrared cutoffs in these continuum 

laws. 

Greek letters a, Aq will be used for various parameters in the models defined in this 

paper. 

Lower case letters will be used for functions which define growth rules for the aggregate. 
These include the functions /, s, t defined latter. 

Latin letters x, z represent points in the complex plane. The number t represents time, 
either as a discrete number of steps or as a real number in a continuum limit. The numbers 9, (p 
represent angles, while the function 9{j) is the step function. 

Unless otherwise specified, subscripts attached to functions will be used to denote deriva- 
tives, thus Fx is the derivative of F with respect to x. As an exception to this, the expression 
f\,e{z) will represent a function parametrized by A and 9. 



2 Conformal Model and Continuum Growth Law 

A model for growth is defined. From this model, the Shraiman-Bensimon equation is derived for 
a function F which maps the unit circle in the complex plane onto the boundary of the growing 
object. Defining 



and making some approximations, we obtain a continuum growth rule, equation (|32D, which still 
includes effects of noise and finite cutoff. 
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2.1 Conformal Model 

The following model for DLA in terms of conformal mappings leads to results that are apparently 
numerically equivalent to DLAQ. It should be noted that the RG in this paper relies upon a 
continuum approximation to this model; this continuum approximation could also have been 
obtained from the lattice version of DLA without reference to the conformal mapping model, 
but the conformal mapping model provides a better justification for the continuum equation. 

In this method, we deal only with the analytic function F, which is defined as the 
analytic function which maps the unit circle in the complex plane onto the boundary of the 
growing cluster. We introduce two parameters, a and Aq, where a = 2 corresponds to DLA 
and Aq is some constant determining the size of an individual random walker. To grow the 
object, first pick a random point x = e*^ on the unit circle. Then, calculate Fz{x), which is the 
derivative of F at this point. We define 

A = Ao(F,(x)F;(x))--/2 (2) 

The case of q = 2 will correspond to DLA and it is that case that will be considered 
from now on; other cases will be briefly discussed in the conclusion. Then in a given growth 
step, F{z) is replaced by 

FifxA^)) (3) 

where / is a function that produces a small bump at angle 6, with linear dimension of the bump 
of the order of the square root of A. X,9 are parameters that define the function /. An explicit 
example of /a,o is given by: 

1/2 



For 6^0, we have fx,e{z) = e'' hfl{e-'' z). 
In the small A case, / reduces to 

z + \z{z + e''^)/{z-e'% (5) 

and by averaging over angle we may determine a continuum growth law. A picture of a cluster 
produced by this growth rule is shown in Fig. 1. A picture of the effect of the mapping / on the 
unit circle is shown in Fig. 2. 



1 + A 

2z 



(z + l)\ z + l + Jz^ + l-2z 



1 - A 
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2.2 Continuum Growth Law 

The continuum growth law for DLA is known to be equivalent to the Hele-Shaw dynamics, which 
obeys the law 

lm{FtF;) = 1 (6) 

where F is a function which maps the unit circle in the complex plane onto the boundary of the 
aggregate and subscripts denote derivates with respect to time or to angle (p on the circle. This 
law may be rewritten as 

ReiFt/zF,) = l/iF^F*) 
Re(Fi^) = 1/|F,| 

where here z = e*"^. Finally, this second growth law may be rewritten as 

F, = F./|(F.,e-»)F;(."r'.^ (8) 

This equation, the Shraiman-Bensimon equation, results from substituting the small A expansion 
of equation (|5|) into the equation (|3|) for the dynamics of F, where F{f{z)) is approximated by 
F{z)+F,{z){f{z)-z). 

In the Shraiman-Bensimon equation, one may divide both sides by zF^ and then take 
the real part of both sides. This will recover equation (|^) and show that the two equations 
are equivalent. Equation ^ implies that the normal velocity of the surface at a given point is 
proportional to the local electric field at that point. 

Equation (|5|) may be rewritten as (taking = for simplicity) 

z + \z{z + l)/{z -l) = z + \z{l + 2/z + 2/z^ + 2/z^...) (9) 

Therefore, the effect of the integration over angle in the continuum growth law is to project out 
negative Fourier components in A considered as a function of angle. The factor of 2 difference 
between the zeroth component and all other components will be important later. 

It will also be useful to define a continuum law for another function G which is defined 
by equation (||). This function has several advantages. The equation for A then becomes 

\ = \oG{x)G*{x) (10) 

which has a simpler form than equation (^) . This has a physical interpretation that G determines 
the strength of the electric field at point x. Also, G is the derivative of the inverse function 
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of F and the inverse function of F has a more natural growth rule than F does. That is, if 
F^'^{F{z)) = z, then under a growth step with given A and x, the function F~^{z) changes into 
/^^(^^^(z)). The continuum law for G, as obtained by using the definition of G and the growth 
law for F, is 

Gt = G,J^Gie'')G*{e'')z{z + e'')/{z-e'') 

-GJ ^G{e^^)G*{ei0)[z{z + e'^)/{z - e% 
It is also useful to consider the continuum growth laws for the power series of F and 
G(z). Writing 

Fiz) = Fiz^ + Foz'^ + F_iz-^ + ... (12) 
G{z) = Goz° + G-iz-^ + G-2Z-^ + ... (13) 
then equation (^) is equivalent to 

{G.j)t = {j -2k-l)Y^ G^kG^iG*_^6{k + l-m- j)2e{j - k) (14) 

k,l,m 

where the discrete step function is defined by 

{1 for j > k 
1/2 for j = k (15) 
for j < k 

An continuum equation may also be written for the power series expansion of F, but we will 
not need to use such an equation. 

There are some problems with directly applying the continuum growth law above, in any 
of its forms, to the discrete random process that defines DLA. The continuum law leads to the 
appearance of cusps in the contour of the cluster after a finite time, and the continuum law is 
deterministic while the discrete law is random. However, the continuum law must have some 
applicability to the discrete cluster growth, because, for example, the conserved quantities of 
the continuum law are approximately conserved by the random process^. Thus, we will try, in 
the rest of this section, to correct the problems in the continuum law so that it may be of some 
use in describing the discrete, random growth process. 
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2.3 Ultraviolet Cutoff 

The above formulation of the problem suggests a method, outlined in this section, of inserting 
an ultraviolet cutoff into the growth law. This cutoff will be inserted by hand, and then the 
parameter of the cutoff will be adjusted to obtain the correct microscopic scale. 

In the discrete conformal mapping model, F never develops cusps because / is always 
well behaved. The specific form of / does not matter; all that matters is that the approximate 
expansion for / given by f{x) = z + \z{l + 2/ z + 2/2;^ + ...) is only correct for small negative 
powers of z. The power series expansion of / is cutoff at some point because / is well behaved. 
This cutoff depends on A, which itself depends on the angle at which growth is taking place. 
The approximation made in inserting the ultraviolet cutoff into the continuum law is that the 
cutoff in the power series expansion for f depends only the average value of A over the circle at 
the time of a given growth step, and not on the exact value of A where the growth is occurring, 
where the dependence on the average value of A is chosen in such a way as to produce the correct 
microscopic scale in the DLA growth process. 

Then a simplification follows. Suppose the regularized form of / is chosen to be 

(l + e)z + e^ 

f{z)x,e = z + z\ ^ (16) 

(1 + e)z — e*" 

where e is a fixed function of (A), which is defined to be the average value of A over the entire 
circle. 

It is worthwhile also to define 

A = l/e (17) 

Then A represents the highest power of z that will occur in the growth law. 

Let us make a change of variable. We will replace F{z) by the function F{{1 + e)^^z) 
and G{z) by G((l + e)~^z). At the same time / is replaced by (1 + e)/((l + e)^'^z). Then the 
continuum law (^) for F becomes 

Ft = F^j ^{F^{{l + e)e^')F:{{l + e)e^')y\'-^ (18) 
where now the cutoff dependence has been moved to the derivatives of F. The continuum law 



(11) for G becomes 

Gt = G,f§^G{{l + e)e'')G*{{l + e)e'')z{z + e'')/{z-e'') 
-GI^G{{1 + e)e'^)G*{{l + e)e^^)[z(z + e^^)/(z - e'^] 



(19) 
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The value of e must now be calculated. 

Before the averaging process, the dependence of the cutoff on A is easy to determine by, 



it is not necessarily the case that the average cutoff e will be determined in the same way from 
the average value of A. The averaging may introduce nontrivial behavior. Instead I will look for 
the dependence of e on the cluster size; since the cluster size and the average of A are related, 
this is an equivalent procedure. 

Expanding the cutoff in the continuum growth law to linear order in e yields an additional 
term in the equation for Ff. This additional term changes equation (^) to 



where z = e^'f' and R is the local radius of curvature. This is a surface tension term. 

The basic idea will be to adjust e to produce the correct size for microscopic features; 
this size is the size of an individual walker in the lattice formulation of DLA. 

A dimensional analysis argument may help understand the size of the cutoff. This dimen- 
sional analysis argument will relate the dependence of cutoff on macroscale to the dependence 
of cutoff on microscale. 

The function F{z) may be assigned the dimension of length, and z may be made di- 
mensionless. This means that we are interpreting 2; as a parametrization of the cluster. Then 
both A and e are dimensionless. We know that A must be a function of the size of an individual 
walker, but then since the dimensional argument implies that A is dimensionless, A must be 
proportional to some power of the ratio of the size of the object to the size of an individual 
walker, as this is the only way to form a dimensionless number. Let tq denote the length scale 
of an individual walker. 

Let us see how to measure the linear size of the object. Recalling the expansion of F in 
Fourier coefficients given by equation (|l^), by a theorem on univalent functions the size of 
the object is at most 4 times Fi, the leading term of the power series. Asymptotically, Fi will 
measure the size of the object. 



for example, expanding the function / as defined by equation (Q). After the averaging process 




(20) 



The additional term may be written as 




(21) 
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To fix the minimum radius of curvature at rg, the size of an individual walker, the cutoff 
e must be chosen so that e/i?, the surface tension term, balances the electric field at the given 
radius. To determine the radius at which they balance, we need to make an assumption about 
the singularities of G. Let us assume G has simple poles and therefore Fz has simple zeroes. 
Suppose we look at points near a zero of F^. Without loss of generality, take this zero to be at 
point zq where zq = 1 — 6 with 5 some small positive number. Locally we find 

Fz^z-zo (22) 

The electric field at z = 1 is proportional to 1/5. The radius of curvature at z = 1 is proportional 
to 1/5^. For the surface tension to balance the electric field we require 

1/S = (23) 

This implies that 5 = e. Then, since R = 1/5^, we need that e = V^. If R is equal to tq, we 
find 

e oc r^^ (24) 

The dimensional argument then implies that e oc (ro/-Fi)^/^, where now the proportion- 
ality constant is dimensionless. This implies that 

A « (Fi/ro)i/2 (25) 

In the actual growth, Fi is changing in time, but rg is constant. Thus the time dependence of 

1/2 

A is determined by A oc . As expected, the cutoff is increasing in time. If the power series 



expansion for G is defined by equation (p~3D, then 

Go = (26) 

so that it is also possible to measure the size of the object using the power series expansion of 
G. 

As a further explanation of the dimensional analysis argument, it may be directly shown 
that, if the cluster is approximately circular, with a small bump on it, then the dependence of 
A on the cluster size is correctly given by equation (^5|). The approximate circularity means 
that instead of simply stating that the electric field is proportional to 1/5 and that the radius 
of curvature is proportional to 1/6'^, we keep track of the proportionality constants in terms of 
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Fi, and then directly show that e oc (ro/Fi)^/^. The advantage of the dimensional argument 
is that it is possible to make this argument without any assumptions on the macroscale of the 
cluster; the dependence of A on rg follows from microscopic considerations, and the dimensional 
analysis argument then yields the dependence on Fi. 

One might worry that for the actual aggregate the poles will not necessarily be simple 
poles. As Halsey et. al. have shown [10|, the surface is described by wedges with a non-zero 
opening angle, and the singularities exist on a fractal set. However, the continuum growth law 
of equation ( [Tl| ) only produces simple poles in G. It is only the dynamics that lead apparently 
to the creation of non-simple poles, via an accumulation of simple poles. Therefore, the cutoff 
will be imposed as if the poles were simple. 

Because the object grows, A increases with time. This is what leads to nontrivial dynam- 
ics and to a fractal dimension less than 2. If instead of varying with time, the cutoff A were held 
constant, then the aggregate would asymptotically grow at a A-dependent rate proportional to 
the square root of time, and would have a fractal dimension of 2. 



2.4 Noise 

From now on, G will be the function of interest and F will be ignored. There are two reasons 
for this: the continuum law for G is simpler, and noise may be more easily inserted into the 
law for G. The actual growth of G is not deterministic; we may write the actual growth of G 
symbolically as follows: actual growth of G = continuum growth of G -|- (actual growth of G 
— continuum growth of G). The term in parenthesis represents noise. This noise term will be 
written as 

sG^ (z) (27) 

where s is a constant with dimensions of inverse time and G^ is some function of z. We 
approximate that G-^ vanishes on average. If we expand G^ in a series as 
we will write the average of G{gJ* as 

{G{Gf;)=6,,N{j) (28) 

where N is some unspecified function. We will also assume that any average of several G^ can 
be written as a product of pairwise averages. These are the essential approximations in the 
noise. 



11 



With noise included, we modify equation (p^) to 

Gt = G./fG((l + 6)e''^)G*((l + 6)e^^)z(z + e^^)/(z-e^^) 

-C/ + e)e^^)G*((l + e)e'^)[z{z + e^^)/(z - e*^]^ + sG^ (z) 

The notation G-'^ is used for the noise because, in the perturbation theory, the function 
G^ will play a role similar to a free field in field theory. 

The approximations are justified for two reasons. Since the noise is essentially generated 
by the dynamics, that is a small amount of noise will be amplified by the continuum growth, 
the dynamics should not be very sensitive to how the noise is inserted. This means we need 
not worry about the exact form of N{j). Second, since the actual growth law for F^^ is rather 
simple, involving a function acting on -F^^, it is most natural to insert the noise into F^^, or 
into G, which is the derivative of F~^. Inserting the noise into the growth law for F, which has 
a more complicated growth law, may have a different effect on the overall dynamics. 

2.5 Continuum Limit In Momentum Space 

In order to make the RG calculations easier, I will also take a continuum limit for the Fourier 
components of G and G-^. This will result in replacing the discrete sums of equation (|T^ ) with 
integrals. This amounts to a change in the geometry of the growth; instead of parametrizing the 
boundary of the growing cluster by a point on the unit circle, we will parametrize it by a point 
on the real line in the complex plane. 

In the neighborhood of a given point on the unit circle, such as z = 1, the unit circle 
is locally approximated by a straight line. As we look at shorter and shorter length scales, this 
approximation becomes more and more accurate. The equation z = e*^ is approximated by 
z = 1 + iO. Thus, on short scales, in the neighborhood of z = 1 we can approximate 

F{z) = Fiz^ + Fqz^ + F^iz-^ + ... « Fi(l + ie) + Fo + J dj F{j)e-'^^ (30) 

G{z) = Gqz^ + G_iz-i + G_22"^ + ...» Go + j dj G{j)e-'^^ (31) 

where F{j) and G{j) should be considered as being defined by the above equations. They are 
defined so that j is always a positive quantity. This approximate form for G{z) will break down 
for 9 of the order of 1 radian. This implies that the Fourier expansions will break down for low 
values of j. This has the effect of an infrared cutoff; the cutoff will be at j of order /i, which 
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is a number of order 1. The cutoff /i is constant in time, unlike A, but under the RG we wih 
find it convenient to rescale /i. After introducing the equation of motion appropriate to this 
approximate expression for G I wiU then explain the effect of nonzero on this equation. 



The symbol j in equations (\3(\), ( 31 ) will be referred to as a momentum since it will play 
a role in the perturbation theory of the next section equivalent to that of a momentum in a 
perturbation theory for a field theory. 

We take equation (^) and transcribe it to this continuum approximation. Noise is added 
as in equation (p9|). The result is: 



GtU) = {l/fi)Jdkt{j,k)G{k)J JdldmG{l)G'^{m)eM-{l + m)/A) ^^^^ 
x5{k + l-m- j)2e{j -k)+ sGf {j) 

where k) is some general function (initially it is proportional to (j — 2k — 1) as in equation 
14), and where the factor of l//x is inserted to produce the correct dimensions for t in the RG, 



as will be clear later. The insertion of the factor of l//i simply amounts to a redefinition of 
t{j,k). The function t will flow under the renormalization group. The exponential term is the 
appropriate version of the ultraviolet cutoff in the continuum limit. We define the continuous 
step function 6 with by the same equation (|T5|). 

The effect of the cutoff /i is twofold: the (5-function in equation ( |32|) has a nonzero width 
of order fi, and hence a finite height. The quantity 6(0) is of order l//x. Also in the definition 
of the noise, equation (|2^) is replaced by: 

{Gf{2)Gf*{j))=6{i-j)N{j) (33) 

where again the (5-function has a nonzero width. 

By rescaling momentum, the cutoffs ^ and A may be changed, but the ratio of the two 
cutoffs will remain constant. The purpose of the RG will be to integrate the upper cutoff from 
A to A — 6A. Then, for the sake of convenience, the upper cutoff will be rescaled back to A. The 
assumption is made that when A is much greater than /x this renormalization does not change 
the essential physics of the system. 



Physically, equation (32) describes the problem of DLA growth in the upper half of the 
complex plane, where the boundary of the growth is parametrized by a point 6 on the real axis. 
The cutoff fj, has the physical interpretation that growth only occurs in a finite width on this 
axis. 
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For use later, let us define a functional W such that 

W[j, Gik),G{l), G* (m)] = (l/fi) J dk t{j, k) G{k) f Jdldm G{l)G*im) exp(-(/ + m)/A) 

x6{k + I — m — j) 

(34) 

Thus, the right hand side of equation @) is W[j,G{k),G{l),G* (m)] + sG^iJ). The functional 
W is linear in each of its last three arguments. 

3 Adiabatic Assumption and Numerical Predictions 
3.1 Adiabatic Assumption 

An important approximation is made, which changes the problem to one of describing an ag- 
gregate which is statistically unchanging in time. In the end, we will describe aggregates whose 
average size and roughness remain constant. 

The cutoff A is slowly changing in time. As the object grows, A changes more and more 
slowly. The exact structure of the cluster at a given time depends upon its growth at all previous 
times, but since the object spends a long time growing with an approximately fixed cutoff, it is 
expected that the structure of the object at a given time t with resulting cutoff A is determined 
only by its growth during previous times t' with resulting cutoffs A' such that A' is very close 
to A. Times t' such that A' is very different from A will be so far in the past that we do not 
expect them to alter the structure of the cluster. 

Furthermore, if the cutoff is fixed, the equations of motion are homogeneous, in the sense 
that up to a rescaling of time and noise, two clusters, which differ only by a change of scale, will 
have exactly the same growth for the same random noise. More precisely, if G{t) is a solution of 
the equation (|29|), with fixed cutoff A and given noise G^ , then, for any number b, the function 
hG{b'^t) is a solution of equation ( p9|) with the same cutoff A and with noise hG^ and s replaced 
by 6^s. 

For a large object it is then reasonable to make the adiabatic assumption that, despite 
the changing cutoff, up to a rescaling of the cluster, the statistical properties of the function G 
at some given time with some given cutoff A are well described by evolving an arbitrary initial 
G for sufficiently long time using the equations of motion with the cutoff held fixed at that value 
A. 
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First, we will analyze the dynamics of G in the fixed cutoff problem, and then we will 
use the adiabatic assumption to relate it to the changing cutoff problem. 

Let us take equation ( |32[ ) and, holding the cutoff fixed, make a dynamical rescaling of 
the function G as it evolves under this equation. After every time step of length dt we will 
rescale G to G(l — sdt). Then, the rescaled G satisfies the equation of motion 

sG{j)+Gt{j) = {l/is)fdkt{j,k)G{k)JJdldmG{l)G*{m)expi-il + m)/A) 

(35) 

x5{k + l-m- j)2e{j -k) + sG-f{j) 

On average, the rescaled G has constant size. 

The amount by which the cluster is rescaled per unit time, s, is the same s referred to in 
the section on Noise. It is simply a matter of notational convenience to chose these two numbers 
s to be the same. Any other choice of s in the section on noise simply amounts to a redefinition 
of Gf. 

The function G before rescaling is growing in time. To determine how rapidly an un- 
rescaled cluster of given size and given, fixed cutoff grows, we may follow this procedure: evolve 
a rescaled cluster using equation (^) with an s chosen such that the rescaled cluster is of the 
desired size. Then, from the value of s needed to maintain the desired size, determine the growth 
rate of the unrescaled cluster. For the unrescaled cluster, the average of = s. 

Under the adiabatic assumption, we can now use the growth rate for the fixed cutoff 
problem, this growth rate being a function of the size of the cluster and the cutoff, and use it to 
determine the growth rate for a cluster with a changing cutoff. To determine the growth rate of 
a cluster of given size in the changing cutoff problem one can determine the cutoff from the size 
of the cluster and then calculate the growth rate of a cluster of the same size in the fixed cutoff 
problem, using the rescaling trick to determine how quickly that cluster grows. 

Under RG scaling, in fact, we will find that the equation of motion changes in such 
a way that s changes; in fact, s may acquire momentum dependence. Before RG scaling s 
will be negative, since F is increasing, causing G to decrease. After RG flow, an appropriate 
combination of t, s, and G-^ goes to a universal value. The adiabatic assumption will mean that 
we assume that at every instant in the original DLA problem, the function G is described by a 
function in the long time limit of the problem with a fixed cutoff and a rescaling term s, where 
s is picked to obtain the correct overall scale for G. 

On average, Gt in equation ( |35|) vanishes. What is left of Gt after the rescaling process 
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is just fluctuations about the average growth. There may in fact be solutions such that Gt 
vanishes identically, but this is unimportant. In the RG, even with Gt non-zero, the s term 
and interaction term {G^ term) will determine the nature of the aggregate. As will be shown, 
under RG flow, the Gt term flows under RG so that, the fluctuations in Gt{j) decrease as j 
decreases. The lower momentum terms then, in the unrescaled problem, will have their growth 
more accurately given just by the sG term. Fluctuations about this overall growth will be less. 



3.2 Numerical Predictions 

It will be worthwhile to mention at this point that already some definite numerical predictions 
can be extracted from the above work. Since the RG that follows relies upon the continuum 
equations, it is good to independently check the validity of these equations for describing DLA. 

If any RG is to hold, the coefficients of G must obey some scaling law. In the last section 
of this paper, such a scaling law will be shown numerically. The coefficients decay with a power 
law. It will be the purpose of the RG to calculate this power law. 

Since the absolute value of G is equal to the local electric field, there is a close connection 



between G and the multifractal exponents of references |12, 13, 14|. These exponents are defined 
by equations (|6S|),([73|). The (2n + l)-th power of the electric field, integrated over the object, is 
given by / d6{G{9)G* (9))"' . This integral over 6 can be converted to an integral over components 
of G in momentum space. Cutting those off at momenta A is equivalent to cutting the real 
space integration off at a length scale of order A^^. Since the multifractal exponents are defined 
in terms of the scaling of powers of the electric field against length, we obtain an equivalent 
definition of multifractal exponents in terms of scaling of powers of G against cutoff. This 
discussion of exponents will be done in more detail later, after the RG permits us to calculate 
these exponents analytically. 

As a check of the adiabatic assumption, a numerical simulation was performed using 
the original discrete conformal mapping model defined in section (1.2). In this simulation, after 
every growth step, the object was shrunk by some constant factor. Visually we could not see any 
difference, in the growing region on the surface, between the cluster shrunk after every growth 
step and another cluster which was not shrunk. A calculation of fractal dimension also failed to 
show any significant differences. 
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4 Perturbation Theory 

A perturbation theory is developed for the equation of motion (|35|). This permits in principle 
the calculation of any correlation function of the theory in terms of noise averages. In practice, 
a resummation of the series is employed which expresses multi-point correlation functions in 
terms of two-point correlation functions. 



4.1 Perturbation Rules 

Using the adiabatic assumption, the equation of motion (|35|) looks very much like the variation of 
an action. Although there no such action can be found, a perturbation theory will be developed, 
based on this analogy, to permit the use of techniques from field theory. 

This technique is very similar to that used, for example, in solving the Navier-Stokes 



equation [15|. Such a perturbation theory has been known for many years. Before proceeding 
with the details, let me summarize the essential attributes: a perturbation theory is developed by 
expanding G in powers of the noise, G-^, and expanding correlation functions of G in terms of two- 
point correlation functions of the noise. Since the noise is amplified by the dynamics of equation 
(|35|) , this expansion is not expected to converge. However, after resummation of the series, it 
becomes possible to replace this by an expansion not in the two-point correlation function of 
the noise, but in the two point correlation function of G itself. Further, a resummation of series 
leads to a resummed propagator (defined below). Unlike the Navier-Stokes perturbation theory 



for turbulence [15 1, it will not not be necessary to define a resummed interaction vertex. The 
above procedure leads to a well defined perturbation series in "skeleton" diagrams. One point 
that will be necessary in the following treatment that is not necessary in the case of turbulence 
is that G is a complex field, and thus the propagators will be directed lines; the notation that 
follows will therefore differ from that seen in Navier-Stokes problems. 

The perturbation theory is constructed as follows: the equation of motion ( [35| ) permits 
us to solve for G(j) in terms of a cubic in G{j). Using the definition of W in equation ( p^ ) we 
write 

G{j) = {s{j) + {w[j,G{k),G{l),G*{m)] + s{j)Gf {j))) (36) 

where s(j) is used instead of s because s may, under RG flow, acquire j dependence. 

The operator (s(j) -|- ■§i)~^ may be expanded as a formal power series in In Navier- 



17 



Stokes perturbation theory, this operator is referred to as the propagator, in analogy with a 
similar object in field theory. 

We may then iteratively solve equation (E^) as follows 



G(i) = ^W[j,G{k),G{l),G*{m)]+GfU)+Oyg-t 



^^W[j,Gf{k),GHl),Gf*im)]+Gf{j) 

^^k, G(n), G(o), G* (p)], G/*(m)] 

+... 

+0(1) 



(37) 



The iterative solution of equation (36), by solving for the values of G{k), G{1), and G*{m) on 
the right hand side, is an expansion in powers of the interaction, t. This is simultaneously an 
expansion in powers of the noise, G^ . The zeroth order and first order terms in t, and one of 
the second order terms, have been written in equation (^). At any point in the process, one 
can stop the perturbation expansion by setting G{j) = G^ (j) plus higher orders in Thus, in 
the perturbation expansion for G one sums at every stage over two possible expressions for G: 

GU) = {l/s{j)){l/f,)Jdktij,k)Gik)JJdldmGil)G*{m)exp{-{l + m)/A) ^^^^ 
x6{k + l-m- j)2e{j -k) + O(^) 

or 

G(i)=G^(i) + 0(|) (39) 

These operations can be represented graphically with Feynman diagrams, in which t is 
an interaction term and s is like a mass term. This leads to a series of diagrams for G. These 
diagrams have no loops. 

The quantities computed as described above depend upon the specific realization of the 
noise. Since we are interested in average quantities, we will compute correlation functions. A 
correlation function is defined as an average over noise of a product of several G{j), G*{j) with 
the same total number of G and G* , and will be written as 

{G{h)Gij2)...G{jn)G*{ki)G*ik2)...G*ikn)) (40) 

When computing these averages, the average over noise must be taken. This is done by taking 
noise terms resulting from the above expansion for G and G* and contracting them with each 
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other in all possible pairwise fashions. Each contraction of two noise terms at momenta j, k 
leads to a factor of 5{j — k)N(j), as given by equation (p3|). This leads to diagrams with loops. 

Fig. 3 indicates how diagrams for the theory are drawn. There are three types of diagrams 
that may be drawn. They are drawn in essentially the same fashion, except that different 
meanings are assigned to the external lines and different numbers of noise contractions are 
included. A diagram for a correlation function has one external line for each term in the 
correlation function to be computed. Diagrams for quantities other than correlation functions 
may have different meanings assigned to the external lines. These diagrams all involve one or 
more external lines being expressed as a function of the other external lines. This may occur 
either as an expansion of G directly in terms of the noise, or as a piece of a diagram that 
occurs inside another diagram. For example, the diagram (a) in Fig. 3 is not itself a correlation 
function, but it represents a term that may occur inside a computation of a correlation function. 
For this diagram, the line for G{j) is expressed in terms of other lines, which in turn may be set 
equal to the noise, or may be further expanded. 

The notation for diagrams is the following: crosses denote contractions on noise, and 
the directions of the lines indicate complex conjugation and orientation within the diagram. A 
line may be said to carry momentum j; when two lines are contracted, they must carry the 
same momentum and all external lines are assigned a momentum determined by the particular 
correlation function to be computed. A t vertex has 2 lines entering and 2 lines leaving, cor- 
responding to G and G* . The G{j) line for a t{j, k) vertex is drawn as entering, the G*{m) is 
drawn as entering since it is complex conjugated, while the other lines are drawn as leaving the 
vertex. The j and k lines will always be drawn parallel to each other in a t vertex; the / and m 
lines will be drawn at an angle. The distinction must be made between j, k and /, m because j 
and k enter into t while I and m are summed over blindly. 

The rules for diagrams for correlation functions may be summarized as follows: to cal- 
culate the average of the product of G{ji)G{j2)---G{jn)G* {ki)G* {k2)-..G* (kn), where ji,ki are 
various numbers: draw one external line for each term in the correlation function. The terms 
in G should be drawn as entering the diagram with momentum ji , while those for G* should be 
drawn as leaving the diagram with momentum ki. Draw all possible diagrams, assigning a factor 
of t{j, k)29{j — k) to each vertex and a factor of l/s{j) + O(^) for each line, while conserving 
the sum of ingoing and outgoing momentum at each vertex. For each noise contraction assign 
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a factor of s'^{j)N{j). Finally, the perturbation theory may give rise to a term such as ^(0). 
This will arise from something like / dj5{j — k)9{j — k). While 6{j) is for negative j and 1 



for positive j, 9{0) will be taken to be equal to 1/2, according to equation (15). This is a result 
of the factor of two difference between powers of z in the expansion for /; this difference was 
discussed in reference to equation (^). Alternately, if we recall that all 6 functions should be 
assigned some nonzero width, then the above integral evaluates to 1/2. 

4.2 Resummation 

This theory exhibits a spontaneous breaking of circular symmetry. One starts the growth process 
with a circularly symmetric cluster, which implies that G{j) = for j > 0, but this is not a 
stable state. Instead, the dynamics evolves G to one of an infinite number of states with nonzero 
G{j), although on average G{j) = for j > 0. If we impose some boundary condition, such as 
G = 1 at time t = 0, and look at G for much later times, the perturbation theory in a small 
noise term causes us to to reach large values of G. Having imposed these boundary conditions, 
the noise is amplified by the dynamics and will grow large. 

This large growth of noise means that if we evaluate the two-point correlation function 
{G{j)G* (k)) we will obtain some answer of the form 6{j — k)\G'^{j)\av, where \G'^{j)\av is an 
appropriate function of j, and \G'^{j)\av is much larger than N{j). In the long time limit, 
\G'^{j)\av is not a function of time. 

Then, we may imagine that, when calculating any other correlation function, at any 
stage in the perturbation theory, two lines which were contracted to obtain the value N{j), can 
instead have their contraction dressed with additional interactions to convert A^(j) to \G'^{j)\av 
Two important points must be made about this procedure. First, it is important that this value 
\G'^{j)\av is completely uncorrelated with any other values of G in the diagram; all correlations 
are already taken into account by the interaction vertices of the diagram. Second, we must 
be careful to avoid overcounting; since each contraction in a diagram includes many diagrams 
involving dressing the contraction in various ways, we must not further dress these contractions. 

For notational convenience, I will continue to write G^ everywhere, but now 

{Gf{j)Gf*ik)) = 6ij-k)\G\j)U (41) 

G^ is similar to a free field in field theory. Any average product of several G^ can be decomposed 
into pairwise products. This is the resummation of noise into two-point correlation functions of 
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G that was mentioned above in the comparison to Navier-Stokes perturbation theory. G-^(j) is 
now an unknown function of j; one of the goals of the RG will be to determine the dependence 
of this function on j. 

Similarly, the function s{j) can be resummed. Any single line between two interaction 
vertices, which would normally be represented by a factor of 1/s, will instead be represented 
by where s^^f takes into account possibilities of dressing that line, without interactions 

with other lines, s^^f is an "effective" s. This is the resummation of propagators. 

As a result of these two resummations, to avoid overcounting we must require that we 
do not count diagrams in which some portion of the diagram that contains interactions has only 
two lines leaving it. Sec Fig. 4 for examples of contributions to j, contributions to s'^^^, and 
diagrams which can not be included in the theory as they would overcount contributions. 

As one will be able to verify after performing the RG calculation, within the formal 
power series expansion of {s + ^)^^- one may neglect all terms in ^ in the RG calculation 
of the next section. Terms only flow to higher powers of ^. Because \G^{j)\av is taken to be 
constant in time, terms with ^ will drop out in many places. 

Then, the perturbation rules, ignoring time derivative, may be summarized as follows: 
draw all diagrams (subject to the rules forbidding overcounting), with appropriate external 
lines, assigning a factor of l/s^^f (j) to each line, and a factor of t{j, k)6{j — k) to each vertex, 
and a factor of {s^^f)'^{j)\G'^{j)\av to each contraction on noise, while conserving momentum 
everywhere. 

4.3 Example Calculation 

It will be useful to give a simple example of applying such a perturbation theory to a non- 
interacting system. For example, consider the following simplified equation of motion: 

Gt{j) + sG{j) = t{j,j)G{j)G{0)G*{0) + sGf{j) (42) 

where t{j,j) is proportional to j — 2j — 1 which is equal to —j — 1. Let us suppose Gf IS 
constant in time, to simplify the problem further. Then, let us define the time scale so that 
t{j,j) = —j — 1. Let us pick the desired scale of the cluster so that G(0) = 1. Physically, G(0) 
is very large compared to the noise G^{0). This mean that for a stationary state (after all, 
s is adjusted to produce a stationary average size) we may let G-'^(O) be small and we need s 
approximately equal to —1. 
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Now that the values of s, t are fixed, we may find the solution, either using perturbation 
theory or using a straightforward solution. The latter method gives —G{j) = —G{j) — jG{j) — 
{j). Therefore, G{j) = G^{j)/j and {G{j)G* {j)) = N{j)/j'^. This noninteracting system is 
stable. 

The perturbation theory for G{0) is slightly tricky. If we ignore the time derivative, we 
can only reach small values of G{0) in the perturbation expansion. However, in fact, for the given 
s and t, a small ^•'^(0) produces a large G(0), under the time evolution. Let us suppose that this 
part of the process has been done, and that we obtain a resummed expression for (G(0)G*(0)). 
This resummed expression is |G^(0)|(j„ = 1. Then, we may ignore the time derivatives and 
obtain the expression for the higher G{j) in a straightforward fashion. 

The perturbation theory expansion to G{j) gives the following infinite sum: 

GU) = Gf{j) + {\G\Q)Ut{j)ls)Gf{3) + {\G\Q)Ut{j)ls?Gf{j) + ... (43) 

This is equal to G^(i)(l - \G'^ {Q)\avt{j) / s)'^ = Gf{j)/j. This perturbation expansion is shown 
in Fig. 5. Each line terminates by setting G equal to the noise. Similarly, the perturbation 
expansion to the correlators gives a product of two infinite sums. This product is represented 
in diagrams by taking the sum in Fig. 5 and noise contracting it with its complex conjugate. 

One may define s^ff{j) for this simple theory; the diagrams of Fig. 5 define the inverse 
of s^^-^ . Therefore 

s-ff(j) = s-tij)\G\0)\a. (44) 
which simplifies equation ( |4^ ) to G{j) = -^^jjjjjsG^ (j). 

5 Renormalization Group Calculation 

We investigate the effect of changing the cutoff in the equation of motion. This leads to the 
introduction of new diagrams to describe the changes in the theory as a result of lowering cutoff. 
It is shown how to incorporate these into a change in s,t. The fixed point is found. 

5.1 Lowest Order Contributions to RG Flow 

In the RG calculation, first the cutoff is lowered from A to A — SA. If the cutoff is imposed in 
a smooth fashion (interaction term t{j,k)G{k)G{l)G* {m) exp[—{l + m)/A]), the change in the 
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theory under a change in the cutoff can be obtained by adding an additional term to the equation 
of motion (^), equivalent to the original cubic term, except that the cutoff exp[— (/ + m)/A] is 
replaced by (5A)f?(exp[— (Z + m)/A])/9A. This new term is 



The sum of the original interaction term plus this new term is equal to the interaction term at 
reduced A. This term will have a circle around the vertex when it appears in a diagram. 

Because this term is small when the external momenta /, m are small, it does not directly 
enter into correlation functions of the low momentum theory. It enters indirectly, in more 
complicated diagrams. We will then consider various such diagrams which include this term, 
and show that, for the simplest such diagrams, these diagrams may be rewritten in terms of a 
change in s and t. 

This procedure has some slight logical differences with other RG procedures. In other 
procedures, the cutoff is often imposed in a sharp fashion at some momentum. Here, the cutoff is 
imposed in a very smooth fashion; I believe this has certain logical advantages. This procedure 
is similar to the technique of counterterms in the original formulation of renormalization in field 
theory. 

First, I will evaluate the simplest diagrams to which this new term gives rise, and do the 
RG calculation, then in the next subsection I will consider other possible diagrams and explain 
why they are neglected. Throughout the RG calculation, when I write |G^(A)|at,, I actually 
mean 



This is just a weighted sum of \G'^{l)\av at / of order A, and due to the smooth cutoff it is this 
weighted sum that will enter into all the diagrams considered. 

The diagrammatic expansion for s*^-^ will now include the diagram shown in Fig. 6. This 



(1///) / dk t{j, k) G{k) J Jdl dmG{l)G*{m){5K/K){-{l + m)/A) exp(-(/ + m)/A) 
x5{k + l-m- j)2e{j -k) 



(45) 




(46) 



changes s' 



■-ff{j) to 




(47) 



The term 29{j — k) in equation ( ^5|) gives 1 in this case, as discussed in the perturbation theory 



rules. 
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Additionally, the new term in the equation of motion can give rise to a diagram as shown 
in Fig. 7a, which can be represented by changing t{j,k). This arises from substituting 

s^ff{l)G{l) = G^{m)t{l, m)G{n)G* {o)5{m + n-l- o)2d{l - m) (48) 

and 

G*{m) = G^*{m) (49) 



into equation (45), taking I of the order of A. The result changes t{j, k) to 

tij, k) + 2i6A/f,)t{j, k)tiA, A + k-j)\GHA + k- j)\ajs'ff (50) 

Assuming A is high momentum and j, k are low momentum, then A-\-k—j = mis high mo- 
mentum. One does not use terms arising from substituting s*^-^ f {l)G{l) = G{n)t{l, n)G^ {m)G*{o), 
instead of the substitution of equation (^8|), because such terms involve too many high momen- 
tum components of G. The diagram corresponding to this term is shown in Fig. 7b. Such 
terms should be ignored, as they will be small when determining the behavior of the system for 
momenta much less than A. Remember that the cutoff is imposed, in the original equation, on 
G{1) and G*{m), not G{k). This means that the presence of the high momentum term G*{m) 
will make such terms small. 

All momenta are now rescaled by j^^g^ to put the ultraviolet cutoff back at A. This 
changes s'^^^{j) to 

s-ffij)-j{ds'ff{j)/dj){6A/A) (51) 

Because of the integration in the interaction term, and the one power of // extracted 
from t{j, k), the dimension of t{j, k) is equal to (momentum)^^. Therefore, under rescaling, the 
function t(j, k) becomes 

t(j, k) - jidtij, k)/dj){5A/A) - kidtU, k)/dk){5A/A) - tU, k)iSA/A) (52) 



Combining the terms resulting from the integration, equation (pOj), with those from the 
momentum rescaling, equation (|52[), every term in the change of t{j,k) either is a function of 
(j — k), or would be a function of (j — k) if t(j, k) were a function of (j — k). For a stationary 
point, then require that 

t{j,k)=t{j -k) (53) 
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If t{j, k) = t{j — k), then similar logic using equations (^7|),(|5l|) requires that s^f^{j) becomes a 
constant. These requirements of momentum independence of s, t likely do not hold at the extreme 
infrared for any true system, or else the equation of motion (equation ( |35| ) using renormalized 
t and s) would have no non-zero solution, but in the scaling region between the infrared and 
ultraviolet, they will hold. 

For the rest of this subsection, the number s is the constant to which s^^ ^ flows under 
the RG, and the number t is t(0). 

At this point, I have done the first part of the RG for the two numbers s and t. It is now 
necessary to rescale G to leave equation (^) invariant under the renormalization and rescaling. 
Let G be rescaled to 

(1 - r((5A/A))G (54) 

This implies a rescaling of s to 

(l + r(M/A))s (55) 

and a rescaling of t to 

(1 + 3r(M/A))t (56) 

Naively one might expect that there could also be an overall rescaling of both s and t by 
a factor a. This would leave G unchanged. Taking derivatives of log(s) and log(t) with respect 
to log(/i) (since the ultraviolet cutoff is lowered to A — (5A and then rescaled back to A, it is 
actually /i that changes in this process), for s and t to be stationary we find: 

r-{K/iJi)t\G\K)U/s + a = {) (57) 

3r + 2(A//i)t|G2(A)|a^/s + a- 1 = (58) 

This in itself does not provide enough information to extract anything useful, because having 
both r and a means that s and t can be scaled to arbitrary values. 



However, in fact, a = 0, as will now be shown. The time derivatives of equation (35) were 
ignored in the RG because they did not effect the flow of s and t. We can still use arguments 
about the time derivative, Gt{j), to show that a = 0. We can use equation (|35| ) to calculate 
Gt{j) of the renormalized, rescaled problem in two ways. One way is to to take Gt{j) for the 
problem before renormalization and rescaling, and then rescale the momentum j and multiply 
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Gt{j) by (1 — r{6A/A)) in analogy with equation (54). Another way is to use the renormahzed, 
rescaled s, t, and in equation ( |35| ) to compute Gt{j). For these two methods of computing 
Gt{j) to give the same result, as required by the fact that they describe the same system, it may 
be shown that 

a = (59) 

A simple way of stating the argument leading to this result is that, although if time derivatives 
are neglected in equation (|35| ) the equation is invariant under multiplying the three terms s, t, 
and G-^ (here referring to the noise before resummation, that is, the noise described by equation 



33) by the same constant, such a multiplication does not leave the time derivatives invariant 
and so should not be allowed in the rescaling part of the RG. 

Then, the fixed point equations (|57|),(|5^) can be solved for 



r = {A/fi)t\G\A)\ay/s = 1/5 (60) 

This result for r is the main result of the first order RG. 

This gives the rescaling of G{j) with j. The magnitude of G{j) must decay as 

Thus 

|G2(i)|a.ocr'/' (61) 

The combination {A/fi)t\G'^{A)\av/s is invariant under a rescaling of G of G by equation (|5^ ) 
and a corresponding rescaling of s and t by equations (|55|) , (|56|) , and gives us the dimensionless 
coupling constant for this problem. From equation (^) the dimensionless coupling constant 
is 1/5, which is not infinitesimal; however, after RG flow, the problem is no longer strongly 
coupled, as the constant is less than 1. 

The reason that a factor of was removed from t{j,k) is now clear; this makes the 
above coupling constant truly dimensionless. As a result of the removal of the factor of 1/;U, 
the dimension of t{j, k) is (momentum) However, \G'^{j)\av has the dimension of momentum; 
this is because, taking G to be dimensionless, equation (|^) gives \G'^{j)\av a dimension inversely 
proportional to the (5-function. The 5-function has dimension of inverse momentum, and thus 
the end result is to make t\G'^{j)\av dimensionless. 

As a further comment on the dimensionality of the coupling constant, recall that the 
(5-functions have a finite height proportional to 1/^. This finite height changes under the RG, 
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which imphes a rescahng of \G'^{j)\av under the RG flow; by multiplying t by we shift this 
rescaling onto t. 

An argument was made above, leading to equation (^), involving Gt for the problem 
before and after RG. We may extend this argument and also say something about the magnitude 
of Gt{j) for different j. The quantitities sG and tG^ remained constant under the RG, as 
a result of the rescalings of G, s, and t and the renormalizations of s and t, from equations 
I) , ([50|) , (^) , (p2|) , (|5^) . We also need Gt to remain constant since this is also a term in equation 
Suppose the characteristic inverse time scale for fluctuations in G{j) is uj{j). Then Gt{j) 
is of order uj{j)G{j). For this combination to remain constant, uj must change as a result of 
the rescaling of j in the RG. In fact, uj{j) must have the same log derivative under RG that s 



does, although the log derivative of s results from renormalization (equation 47) while the log 
derivative of a;(j) results from rescaling. This implies that 

a;(j)oc//5 (62) 

This means that for smaller j, the time scale for fluctuations is longer. Returning to the original 



problem, as described by equation (29) with a time-dependent cutoff, this means that the lower 
Fourier coefficients grow at a roughly constant rate. This self-consistently justifies the adiabatic 
assumption of section 3. 



5.2 Other Contributions to RG Flow 

One may imagine that the new term of equation (^5[), representing the rescaling of the cutoff, 
may enter into additional diagrammatic contributions. Various possibilities are shown in Fig. 8. 
I will show that, for low momentum behavior, these terms are unimportant and then discuss in 
more generality why other contributions are negligible. 

One may check by hand that the first example is small if external momenta are much 
less than A. The second example will be discussed below. The third example vanishes due to 
phase space factors. The fourth example vanishes due to phase space factors if the two lines 
leaving the top of the diagram are close in momentum; this means it does not alter the RG flow 
of t{j, k) when j = k. In the rest of the section, various other diagrams will also be said to 
"vanish"; this will only mean that they vanish when considered either at low momentum or, if 
they contribute to the RG flow of t(j, k), when considered at j = k. 
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The fifth example, a contribution to the six-point function, will be seen below to be small 
when calculating correlations of only 4, and not 6 or more, G{j). The sixth example should 
not be considered when the rescaling of \G'^{j)\av is taken into account; since (A//i)|G^(A)|a„ is 
stationary under RG flow, such a diagram is canceled by the various rescalings. 

In order to indicate in general why such contributions may be neglected, I would like to 
define some additional terminology to describe certain paths and contractions in these diagrams. 
When considering a contribution to s*^^ ^ , one may follow one line through the diagram as follows: 
start with the incoming line. At every t vertex, if one enters with G'(j), follow out along G{k), 
not G{1) or G*{m), where the roles of j, k, I, m are as in equation (|35D. If the line one is following 
is contracted (this will be referred to as an exceptional contraction) with a G*{m) leaving a t 
vertex, follow out along G{1) of that t vertex. This path will be referred to as the main line. 
Now, any diagram that includes an exceptional contraction, such as the diagram of Fig. 8a, will 
be small for low momenta, since the G{h) leaving such a diagram will only have a small range 
over which to integrate. For Fig. 8a the main line is given by following the horizontal arrow 
along the bottom of the diagram from left to right, through the t vertex, until it bends up and 
left. Then go down and left through the noise contraction into the t vertex, and then leave the t 
vertex along the line going up and left. Follow this line through its bend back to the right until 
it leaves the diagram. 

The smallness, of the contribution to s^^^ given in Fig. 8a, for small momenta does not 
completely justify the neglect of such terms. For example, when evaluating the RG contribution 
to t, the value of s^^^ used is s^^f{K), not the low momentum s^^^, and thus a high momentum 
contribution to s^^^ may change the low momentum renormalized t. However, even for a calcu- 
lation of s'^^f (A) , the exceptional contraction will mean some reduction in available phase space 
over which to integrate. 

One may check that contributions to s^f^{j) like the third example in Fig. 8 will always 
vanish, regardless of what j is, due to the 6 functions in equations (|35|) , (|45|) . The lines coming 
off of t vertices connected to the main line must be contracted within themselves, not between 
different vertices. In this diagram, the main line is simply the entire horizontal arrow running 
along the bottom of the diagram. Thus, the only contribution to s'^^^ will be the contribution 
of Fig. 6. 

For contributions to t(0), we may define 2 main lines. One can follow the main line of 
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the G{1) or the main hne of the G*{m). These are the hnes one follows if one starts on the line 
for G{1) or for G*{m) and follows through the various contractions as defined above. These two 
lines join at some point in a noise contraction. 

The second example of Fig. 8 has an exceptional contraction and may be ignored. The 
main line starting with the G{1) line leaving the circled t vertex proceeds up and left, then turns 
right, going straight across the diagram to the end. Then it turn down and left, up and left 
through an exceptional contraction into a t vertex. Then it goes down and left into a noise 
contraction, where it ends. The main line starting with G*{m) starts at the circled t vertex and 
proceeds up and right until it terminates at the noise contraction. Fig. 8b is very like Fig. 7b, 
except an additional t vertex has been added to the diagram. If this vertex were removed, this 
diagram would be small for small external momenta. With the t vertex on the diagram, the 
diagram is very difficult to evaluate since to evaluate it requires a knowledge of all t(j, k), not 
just t(A, A). However, the diagram is not only next order in the coupling constant, but also 
small due to the various exponentials present, as may be verified. Any diagram for t with such 
an exceptional contraction will have the same problems. That is why we will ignore them. 

If the main line of the G* (m) has lines leaving it which contract against lines leaving the 
main line of G{1) then the diagram will again vanish due to the 9 functions. A contribution to t 
cannot have both external lines leaving the main line of G* (m) . Therefore both lines must leave 
the main line of G{1), and the line G*{m) cannot be dressed by any interaction vertices. Again 
due to ^-functions, the main line of G{1) can only include one interaction vertex which has both 
external lines on it. For example, the diagram of Fig. 8d has external lines coming off different 
interaction vertices and vanishes if the two lines are close in momentum. 

All that remains is to justify the neglect of six-, and higher-,point functions, such as 
could appear from the diagram in Fig. 8e. If we wish to calculate a correlation function of four 
G{j), and somewhere in one diagrammatic calculation we have a six-point function, some of the 
lines leaving the six-point function must be contracted against each other. Then, somewhere 
in the diagram for the correlation is a contribution to the four-point function which includes 
the six-point function within it. Therefore, the renormalization procedure would have yielded 
this contribution to the four-point function as a change in t. But, we have already, as outlined 
above, obtained all the contributions to the change in t. Therefore, there is no such diagram. 

Finally, the lowest order contributions considered in the previous subsection have a 
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certain universality; considering only the diagrams of that subsection, the nature of the fixed 
point does not depend on the initial form of t{j, k). Higher order RG contributions will depend 
on the initial form of t{j, k). 



6 Fractal Dimension and Multifractal Exponents 

It is now possible to begin extracting exponents of the original DLA model. Different exponents 
correspond to different correlation functions of this model; it will be the purpose of this section 
to determine how to compute exponents from correlation functions. This process depends on the 
discussion of the adiabatic assumption and the assumption used to introduce the cutoff into the 
continuum equation. From those assumptions an unambiguous means of determining exponents 
from correlation functions is given. 

In any actual simulation, there is an ultraviolet cutoff A determined by the ratio of 
macroscale to microscale. In the RG, a power law decay was found for G{j). Since the RG 
describes a fixed point in the scaling region, within the RG itself this power law holds for 
arbitrarily large j. Within an actual simulation this power law will fail at j of order A, where A 
is the cutoff of equation ( |25| ) resulting from the finite size of the cluster in the simulation. Thus, 
in the calculation of exponents that follows, although all correlation functions are calculated 
using the rules of the RG and of perturbation theory, the integrals over correlations functions 
that we will use must be cutoff at momenta of order A, as will be done. 

Since the size of the object follows a power law behavior given by 

Fi oc t^/^ (63) 

where D is the radius of gyration dimensions of the object, we have 

dlog(Fi) 



dt 



oc l/t (64) 



As a side point, strictly speaking this requires that dFi/dt can be replaced by the derivative of 
the average value of Fi, but both numerical evidence and the RG flow of frequency resulting 
from equation (|6^) justify this assumption. However, 

dFi/dt = {X)Fi (65) 
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where (A) is defined to be the average value of A over the unit circle at a given time. Equation 
( pBD may be derived by using equation ( p^ to calculate dGo/dt and then equation (^) to relate 
this to dFi/dt. Combining equations (p^) , (|6^) we get 

(A) oc l/t (66) 

Equation ( p^ ) is equivalent to the electrostatic scaling law first derived by HalseyQ. In the 
continuum approximation, (A) = / dj{G{j)G* (j)) . Here it must be understood that while in 
the perturbation theory this expression is formally infinite, since {G{j)G* (j)) = 6{0)\G'^{j)\av, 
in the above average we remove this factor of (5(0)[17|. Calling D the fractal dimension, 

Ff ^ oc l/t oc (A) oc Gl 1^ ^ ^-2^3/5 ^ ^-2^3/10 (g^) 

This gives the result that 

D = 2 - 1/2 + 1/5 = 1.7 (68) 

The first proportionality in equation (|6^) followed from the radius of gyration definition 
of the dimension. The second followed from the electrostatic scaling law. The third followed 
from the expression for (A) in terms of G and from the scaling of G derived in the RG. The 



fourth followed from equation (26) and from doing the integral. The fifth followed from the 
functional dependence of A on Fi as given by equation (^). 

This is the simplest way to derive the fractal dimension from the above work. The calcu- 
lation of the growth rate from equation ( |66D is essentially a determination of the unrenormalized, 
unrescaled s in equation (|35D. It may also be possible to repeat the same result by using the 
rescaling of s under the RG to obtain the rescaling of the growth rate under a shift in A. 

The multifractal exponents T{q) are defined by 

r(g)=lim(^log(^F5(i))/log(/)^ (69) 

where the surface of the cluster is covered with intervals of length /, and E'^{i) is the integral 
along the ith interval of the gth power of the electric field. Numerical calculations of these 



exponents can be found in references 18]. 



One can try to compute higher multifractal exponents using the RG (the work above 
amounts to computing r(3) and showing that D = t(3)). For example, the scaling of t(5) can 
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be determined by calculating the scaling of 

r-A 



djdkdldm{G{j)G{k)G*{l)G*{m))6{j + k - I - m) (70) 

against the upper cutoff A. This is because the given integral is equal to the desired power of 
the field integrated over the surface of the object. If all the terms in this integral contributed 
with the same phase, the integral would scale as A~^/^~^^. Of course, the terms are independent 
and this misestimates the exponent. It is necessary to use the perturbation theory to evaluate 
the 4 point correlation function. The simplest possibility is to use Gf as an estimate for all the 
G in equation (70). The diagram for this is shown in Fig. 9a. The only terms that would then 



contribute would be when a G and a G* were at the same momentum and the integral would 
scale as A"^/^"^^. The different scaling results from having a different number of momenta to 
integrate over. Another possibility (this is analogous to a tree approximation for a scattering 
problem) is to substitute for the highest momentum G in terms of a f vertex, leaving a six-point 
correlation function, and then take all six G, G* to be Gf. The scaling is then as ^-^1^+^. The 
diagram for this is shown in Fig. 9b. Since this scales more strongly with A, it will be dominant 
in the limit needed to compute T{q). 

The following is the rule for calculating multifractal exponents: Lei n he a positive integer. 
Calculate the integral over 2n-point correlation functions defined by 

dji dj2... djn dki dk2... dkn ^^^^ 

{G{h)G{j2)-G{jn)G*{ki)G*{k2)...G*{kn))6{h + J2 + ... + jn " ^1 " ^2 - ... - kn) 

If this integral behaves, in the limit of h. taken to infinity, as A", where a is some number, then 
r(2n + 1) = (2n) — a/2. The factor of (2n) is the trivial scaling that would result even for a 
nonfractal object; the factor of a/2 results from the dependence on A and from the square root 
dependence of A on length scale. 

In general, we can always find, for r(g), a tree diagram that scales like A~'^^'^~'^^/^'^^'^^'^\ 

Then, 

r{q) = (<? - 1) - l/2(-2(g - 2)/5 + {q - 2)) (72) 
Alternately, another definition of exponents is 

Dg = T{q)/{q - 1) (73) 
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Then 



which is equivalent to 

0.7g-0.4 



= •„ / (75) 



7 Comparison With Numerics and Discussion 

The theory is compared with numerics, and further tests of the theory are proposed. 
7.1 Comparison With Numerics 

In previous work we found that the alternate formulation of DLA using analytic functions |^] 
produces clusters with appearance and dimension similar to those of clusters grown using the 
lattice formulation of DLA. As far as we can tell, the two formulations are equivalent when 
a = 2. 

The simplest comparison with numerics is the dimension itself. 1.7 is very close to the 
accepted value of 1.71. 

Equation (^) for higher multifractal exponents is the same as the formula obtained with 



a wedge model by Halsey et. al.[|lO[, except that the wedge model left the quantities 0.7 and 0.4 
as unknown constants to be fitted to numerics. They define quantities / and a, the dimension of 
the set on which the wedges exist, and the strength of the singularity (hopefully, the reuse of the 
symbol a will cause no confusion), and show that Dg = ■ ^ numerical fit gave a = 0.705, 
/ = 0.42, while a comparison with equation ( [75| ) gives a = 0.7, / = 0.4. It is now known that 
such a simple scaling law is not valid for large <?[|l8|, and in the original paper of Halsey et. al. 
it was suspected that such a law would not hold. 

The possible difference between theory and experiment here for large q should not be 
construed as a flaw in the presented work. First, the above calculation is only a lowest order 
calculation. To higher orders, we may find a form for t(j, k) which has nontrivial behavior. This 
may alter the results from the tree approximation to the correlation function used to compute 
the exponents. Second, we may find interesting behavior if we consider other contributions to 
the correlation exponents, beyond the tree diagrams used above. Third, although the neglect 
of the appearance of 6-point function was valid when considering the renormalization flow of 
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s and t, as discussed in reference to Fig. 8e,, such a neglect is not valid if one actually wishes 
to compute 6- and higher point correlation functions. Fortunately, such multi-point interaction 
terms are captive variables, in the sense that if one knows the behavior of s, t under RG flow one 
may systematically determine the higher interaction terms that will appear. Fourth, the above 
derivation of multifractal exponents involved expressing the exponent in terms of correlation 
functions; this is only possible for odd multifractal exponents. Thus, in fact it is not possible to 
say anything about even exponents in any simple fashion. 

Additionally, there exist some difficulties in numerical calculation of higher multifractal 



exponents. According to the branched growth theory of DLA |16], the time required to compute 
higher exponents is superexponential in the order of the desired exponent. Thus, the exact 
values of the larger exponents may not be precisely given by the numerical experiments. This 
mathematical difficulty may be the source of the controversy which appears to exist between the 
different numerical calculations of these exponents. For example, the value quoted for r(3) by 
Ball et. al. is less than 1.6, which is definitely at odds with the electrostatic scaling law (believed 
to be exact from various numerical calculations), and with other numerical calculations of this 
exponent. The electrostatic scaling law says that t(3) = D = 1.71. There also exists controversy 
about the precise value of the dimension of DLA, as mentioned in Ref. Thus, in fact, it 
is not clear exactly how large is the discrepancy between the above results and the numerical 
results. 

It is also of interest to numerically check the scaling of G{j). This was check for two 
cases. First, for a single cluster as described in the next paragraph; second, for an ensemble of 
clusters as described in the paragraph after that. 

When / dj G{j)G* {j)e~^^^ is plotted against A*^'^, where in reality the integral is a 
discrete sum, one expects to find a straight line behavior. This is what is found, as shown 
in Fig. 10, except that for large A the curve flattens out, since the numerical calculation only 
included a finite number of terms. Also for small A, the curve flattens out at F^'^, which, in 
the long time limit, is vanishingly small compared to the full integral. For the finite cluster 
size of our simulation, F^'^ is not negligible. The clusters here were grown using the conformal 
mapping technique outlined previously. The coefficients of G were computed with a numerical 
Fourier transform, by mapping a large number of points on the surface of the circle (in fact, 
slightly outside the circle, to improve numerical behavior) to the surface of the aggregate (again. 
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to slightly outside the surface of the aggregate) and analytically calculating the derivative of the 
mapping for each point. This technique is not very efficient for growing large clusters, at least as 
presently implemented. It requires 0(A/"^) time to compute N growth steps, but it is very easy 
to calculate coefficients of G using this program. I only used aggregates of around 7000-10000 
walkers. 

As another check, 50 clusters of 6000 steps were simulated, and for each cluster, the 
coefficients of G were computed. The squares of these coefficients were scaled by G(0), the 
overall inverse cluster size, and then averaged together. In Fig. 11 the mean squares of G{j) 
are plotted against j on a log- log plot. Numerical difficulties made it impossible to accurately 
extract the slope in the scaling region. This scaling region extends from j = 5 to j = 35, or 
from log(j) = 1.6 to log(j) = 3.5. Theory predicts that this slope is 2/5 = 0.4. The numerical 
slope is between 0.3 and 0.5, using a least squares fit. The theory line is drawn onto the graph. 
As an additional check, another ensemble of clusters was simulated, with a different Aq and a 
different number of steps. Within the scaling region of that simulation, the slope of the mean 
squares behaved in the same fashion, and, additionally, the mean-square of G{j), after scaling 
by cluster size, for given j, was the statistically the same for the two simulations. 

7.2 Discussion 

A theory has been presented based on the conformal nature of various Laplacian growth pro- 
cesses. A series of approximations were made that produced a modified continuum equation 
of motion; it is hoped that such an equation describes DLA, but even if it does not, it does 
describe some form of nontrivial Laplacian growth. A perturbation theory was developed for 
this equation, and rcsummed. To determine various terms in the perturbation theory, it was 
then necessary to use a renormalization group calculation. This has only been carried out to 
lowest order. It is a peculiar feature of this method that next order calculations are vastly more 
difficult than lowest order calculations, thus as yet there is no analytic calculation of higher 
order effects. Finally, the assumptions leading to the modified model were reversed, leading to 
calculations of quantities for DLA. 

It would be worthwhile to look more closely at higher order corrections, if not analytically, 
at least qualitatively, to see what may happen. To lowest order, t{j — k) flows to an everywhere 
positive function. Using the lowest order t to compute the effect of higher order corrections will 
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tend to lower the value for the dimension predicted by this theory. However, it is possible that 
in a more careful next order calculation, the interaction t{j — k) flows to a function which is 
negative for large (j — k), possibly increasing the predicted dimension. As mentioned above, 
next order effects depend in some way upon the initial functional form of t{j, k). Lowest order 
effects do not. 

Unfortunately, it is not possible to carry out a stability analysis of the fixed point of the 
lowest order RG. All that may be said from the above calculation is that if a fixed point exists, 
other than a trivial fixed point for which t goes to zero, then this fixed point is described by 
this RG. 

It would also be worthwhile to try to extend this technique to other Laplacian growth 



models, such as the dielectric breakdown model. For the dielectric breakdown model|ll|, dif- 
ferent values of r] correspond, in the continuum limit, to different values of a in the conformal 
mapping model of section 2.1. The difference between the dielectric breakdown model and our 
model is that, away from the DLA case, our model uses the same growth probability over the 
surface and varying walker size, while the dielectric breakdown uses a varying growth probability 
and constant walker size. Although this alters the scale of the cluster in our case, we would 
expect the fractal dimension of the cluster grown at a given a with the conformal model to be 
the same as the dimension obtained from the dielectric breakdown model with r] equal to a — 1. 

One might naively try to apply the technique above to the case of a different from 2, by 
replacing equation (|l|) with the definition G = Fz "^^"^ , and using an equation of motion similar 
to (^), with different initial t[j, k). This would lead to physically absurd results, and is in fact 
different from defining G = F^^ and using a modified equation of motion as described in the 
next paragraph. The difference is in how the noise term is inserted. It is important for the 
perturbation theory that products of G^ may be pairwise decomposed, and this property means 
different things depending on whether G = F~^ or G = Fz In the stochastic problem, each 
growth step produces a simple pole of F^^ inside the unit circle; the angular coordinate of the 
pole is random, the radial coordinate is determined by A. In the continuum limit of section 
2.5, the angular coordinate becomes the real value of 9, while the radial coordinate becomes 
the imaginary value of 6; the interior of the unit circle is replaced by the lower half plane. One 
may show, using Cauchy's theorem, that randomly inserting simple poles produces a pairwise 
decomposition property for the random noise in F~^. Therefore, the equation (||) is the best 
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definition of G. 

One can handle the case of a different from 2 by using a modified equation of motion, 
although this may be difficult if a is not even. One would modify equation ( ^5| ) by including 
higher powers of G in the interaction term. For example, for a = 4, the interaction term would 
be of the form 

/ dk t{j, k)G{k) JJdldm dn do G{l)G* {m)G{n)G* (o) 

(76) 

exp(— (/ + m + n + o)/ h)5{k + l + n — m — j — o)26{j — k) 

This would probably be the most worthwhile test of the calculations of this paper; although 
the calculation for a = 4 is far more difficult than that for a = 2, it may still be tractable to 
lowest order. It would not be appropriate to attempt such a calculation in this paper. A few 
preliminary calculations show that one obtains at least the physically correct result that the 
dimension of the a = 4 model is less than that of the a = 2 model; as yet the exact value at 
a = 4 is not calculated [p^ ]. 
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